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Abstract 

The OPERA collabotation has reported evidence of superluminal neutrinos with a 
mean energy 17.5 GeV ranging up to 50 GeV. However, the superluminal interpreta- 
tion of the OPERA results has been recently refuted theoretically by Cherenkov-like 
radiation. We discuss a loophole of this argument from the kinematical viewpoint and 
find it possible to avoid the Cherenkov-like radiation of the OPERA neutrinos. The 
key idea of our argument is to admit the fact that the neutrinos travel faster than the 
observed speed of light while they do slower than the true speed of light in vacuum 
so strictly speaking they are not superluminal but subluminal. Moreover, we present 
a model where these two velocities of light can be constructed by taking account of 
influences from dark matters near the earth. 



1 E-mail address: ioda@phys.u-ryukyu.ac.jp 



1 Introduction 



The OPERA experiment has recently announced a remarkable discovery that neutrinos from 
CERN to Gran Sasso travel faster than the speed of light [1]. Soon after this sensational 
announcement, there have been a lot of activities to attempt to understand this interesting 
phenomenon of superluminal neutrions from the theoretical viewpoint. As a consequence 
of such theoretical studies, it has turned out that there are two theoretical challenges to 
the OPERA experimental results, which are the bremsstrahlung effects [2] and pion decay 
[3, 4, 5]. Since the two problems stem from the same cause that the OPERA neutrinos are 
superluminal, we focus on only the bremsstrahlung effects in this article. We will see that the 
problem of pion decay is also automatically solved by our theory. 

Let us start with an explanation of what we call, the bremsstrahlung effects or Cherenkov- 
like radiation of superluminal neutrinos. On the way from CERN to Gran Sasso, the very 
effect of superluminal propagation of neutrinos would have caused some distortions in the 
beam of neutrinos owing to Cherenkov-like radiation and severely depleted the higher-energy 
neutrinos, thereby making it impossible to observe neutrinos with more than 12.5 GeV energy. 
This theoretical result is obviously against the OPERA results where a lot of high-energy 
neutrinos above 12.5 GeV are observed as seen in the OPERA paper [1] "Data were then 
split into two bins of nearly equal statistics, including events of energy lower or higher than 
20 GeV. The mean energies of the two samples are 13.9 and 42.9 GeV." 

More recently, stimulated with the above theoretical observation, ICARUS group has 
analyzed their data and found that "the neutrino energy distribution of the ICARUS events in 
IAr agrees with the expectations from the Monte Cairo predictions from an unaffected energy 
distribution of beam from CERN. Our results therefore refute a superluminal interpretation 
of the OPERA result according to the Cohen and Glashow prediction for a weak current 
analog to Cherenkov radiation" [6] . 

In order to make a physically plausible theory of the OPERA results, one has to overcome 
this issue at all events. It is our aim in this article to present one solution to this problem. To 
do so, it is worth keeping in mind that since physics of Cherenkov-like radiation of superlu- 
minal neutrinos is purely in the kinematical regime, one has to propose a not dynamical but 
kinematical solution to this problem which would forbid this process from the kinematical 
reasons. 

In addition to it, we cannot help accepting the fact that according to Cohen and Glashow 
if a neutrino could travel faster than the speed of light in vacuum as claimed by OPERA, 
superluminal neutrinos rapidly lose energy via pair bremsstrahlung so few neutrinos with 
energies in excess of 12 GeV are detected at Gran Sasso. 

In this article, let us suppose that the OPERA results are correct at any rate even if 
further experimental scrutiny is surely needed. Then, a question naturally arises: "Is there 
a way out to reconcile the OPERA results which insist on the superluminal propagation of 
neutrinos, with the bremsstrahlung effects which refute a superluminal interpretation of the 
OPERA results by the kinematics?" Below we wish to present a solution to this question. 
The key idea is to introduce two kinds of speeds of light in such a way that the OPERA 
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neutrinos are not superluminal but subluminal under the true velocity of light. 

This article is organized as follows: In the next section, we will review the work by Cohen 
and Glashow. In Section 3, we explicitly make a model which has two speeds of light because 
of interaction between the gauge field and a Galileon-type of scalar field. Section 4 is devoted 
to discussion. 



2 Review of bremsstrahlung effects and its resolution 

Let us begin with a review of the bremsstrahlung effects of superluminal neutrinos [2]. This 
study is based on an old analysis of Ref. [7] where many effects of violation of Lorentz 
invariance have been explored in the framework of Standard Model. 

The basic idea behind the work by Cohen and Glashow is that a neutrino traveling faster 
than the speed of light loses energy by emitting something owing to weak interaction although 
neutrions do not carry electric charges (So this process is dubbed Cherenkov-/i£;e-radiation as 
well). The most dominant decay process of a muon neutrino, which mainly constitutes the 
OPERA beam, is found to be — > + e + + e~ whose threshold energy reads 

atyi 

E = 6 « UOMeV, (1) 
iv? — 1 



for the OPERA neutrinos. 

In the high energy limit where mass of an electron and neutrino can be neglected, the rate 
of bremsstrahlung pair emission T and its associated rate of losing energy are calculated 
to be 

r = k '-^-E\vl-i)\ 

192^ 3 v " ; ' 

^ = - k J3_ E %vl ~ I) 3 , (2) 
dx 192tt 3 v " ; ' v ; 

where k = ^$,k' = ^ are numerical constants. Then, the mean fractional loss of energy at 
each step of the bremsstrahlung pair emission is given by 

-— k 

dx - — w 0.78, (3) 



r# k' 

which implies that about | of the neutrino energy is reduced at one step of emission. 

Moreover, by integrating ^ over a distance L, the final neutrino energy E is expressed in 
terms of its initial one E by 

jL = .1 + hk-^-Avl - ifL. (4) 
E 5 El 192tt 3V ; v ; 

By applying this formula to the case of the OPERA experiment, the terminal energy of a 
neutrino at the OPERA detector is about 12.5 GeV. In other words, few neutrinos reach the 
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detector with energies in excess of 12.5 GeV. Unfortunately, the OPERA detector observes 
neutrinos with the mean energy of 17.5 GeV ranging up to 50GeV [1], so the calculation rules 
out an interpretation that the OPERA neutrinos are super luminal. Incidentally, the ICARUS 
group has recently reanalyzed their data obtained in 2010 and found that no Cherenkov-like 
radiation has been detected in ICARUS [6]. 

Where is there a loophole in the above argument? It is true that we cannot help accepting 
the conclusion by Cohen and Glashow as long as a neutrino is superluminal. Thus, the only 
way out is to consider that the OPERA neutrinos are not superluminal but subluminal, but 
then how can we make them stay subluminal, which appears to be against the OPERA report? 

As an implicit assumption of the argument by Cohen and Glashow, the speed of light is 
considered as a universal constant and does not receive any influences from the environment. 
However, it has been already pointed out that there is the possibility that the velocity of 
the photons might be dependent on energy from the HERA Compton polarization data [8]. 

2 It might be therefore reasonable to conjecture that the speed of light would receive some 
influences from surroundings and consequently the observed speed of light in vacuum might 
become smaller than the true speed of light in vacuum. 

Put differently, if the group velocity v u of the OPERA neutrions is between the observed 
velocity Co of light and the true one c, i.e., Co < v v < c, the OPERA neutrions can be regarded 
as subluminal since the property of superluminality or subluminality is defined by using the 
true velocity c of light in vacuum. 

Then, the problem of the bremsstrahlung effects is converted to a different problem: Can 
we construct a model which has two kinds of velocities of light, which satisfy the relation 
c < c, without conflicting with special relativity? In the next section, we shall present such 
a model based on the recent works [9, 10, 11, 12, 13, 14]. 

3 A model of two velocities of light 

We shall begin with the Lagrangian density of our theory : 3 

where M* is a mass scale which controls the strength of the coupling between the Galileon- 
type of scalar field n [15] and the abelian gauge field A^. The gauge field strength F^ v is 
defined by F^ u = d fM A u — d u A fl as usual. In the absence of the 7r field, the photon travels at the 
speed of light in vacuum since the gauge field satisfies the conventional Maxwell equations. 
In the theory at hand, it is essential to define this velocity, which we denotes ctS C« clS the true 
speed of light in vacuum, which appears in various formulae in special relativity. Thus, the 
causal structure should be determined on the basis of this velocity c and the flat metric r}^. 

2 We identify the velocity of the photons with that of light. 

3 We make use of the flat metric jj m „ = diag(—l, +1, +1, +1) for raising or lowering indices. 
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Now it is easy to rewrite (5) as 



L = ~\W a ~ j^d^n)^ - -^d^d^F^, (6) 



thereby making it possible to read out an effective metric 

f 



ffil) = rT - (7) 



along which the photon propagates. 

Next, let us consider the spherical symmetric configuration for the scalar 

where a is a constant [11]. Note that this configuration breaks the Lorentz invariance spon- 
taneously. Then, it turns out that the effective space-time on which the photon propagates 
has the line element 

ds 2 = g^^dx^dx" = -dt 2 H ^-^—^dr 2 + r 2 dVt 2 2 , (9) 

where dVi\ = d6 2 + sin 2 9d<p 2 . 

It is then straightforward to derive an effective velocity c(r) of the photon for the fixed 
angles (9,<f>) 

c(r)^=(l-4^l% (10) 
K J dt \ M$r 4 J K J 

where for conveniece the true velocity c of light in vacuum is recovered. Note that at the 
spatial infinity r — >■ oo, the effective velocity coincides with the true velocity c, that is, 

lim c(r) = c. (11) 

Furthermore, note that the difference between c and c(r) is positive-definite as long as the 
constant a is non-zero 

a 2 1 

c-c(r) = — — — A > 0, (12) 

which holds when ^<1. 

On our earth, the observed velocity of light is given by c(r) when the photons are located 
at the place whose distance in the radial direction is r from the center of the earth. If we 
assume that the speed of a neutrino, denoted as v v , takes a definite value and is smaller than 
c but larger than c(r) 

c(r) < v u < c, (13) 
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we would have a superluminal neutrino for the velocity c(r) as observed in the OPERA 
experiment 



P(r) = ^— ^ > 0, (14) 
c[r) 

whereas we have a subluminal neutrino for the true velocity c 

^ = !^ <0 . (15) 

c 

Recalling that the property of superluminality or subluminality of neutrinos is now defined 
by using the true velocity c of light in vacuum, the OPERA neutrinos are actually not 
superluminal but subluminal! Hence, we do not have Cherenkov-like radiation for the OPERA 
neutrinos at all since we can always take the rest frame for the subluminal neutrino. This is 
our resolution to the problem of the bremsstrahlung effects of the OPERA neutrinos. It is 
worthwhile to notice that our solution is purely kinematical as desired. 4 

Let us show that the results of OPERA and SN1987A give us information on the mass 
scale M*. First, let us note that the OPERA result yields a condition for the dimensionless 
quantity (5{r) 

P(Re) = V ^ C{R ^ « 2.5 x 1(T 5 , (16) 

where is the radius of the earth and takes the value = 6.4 x 10 8 cm. 

Next, let us utilize the fact that neutrinos from SN1987A to the earth travels at almost 
the same velocity as the true velocity of light in vacuum, so we have a relation 

v v « c. (17) 

With the help of this relation (17), Eq. (16) can be cast to the form 

x ik-k- (18) 

At this stage, we may assume the scalar field tt is sourced by the trace part of the energy- 
momentum tensor on the earth which does not include the contribution from the scalar field 
7r itself. With this assumption, the numerical constant takes the form in the static case [11] 

a=^ = M e L„ (19) 



where M e is the mass of the earth and we have defined the length scale by L* = -77 



4 Rccently, a similar resolution has been also proposed by Nakanishi in Rcf. [16]. It is a pity that his 
interesting article is written only in Japanese. 
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Using Eq's. (16), (18) and (19), is calculated as 



( 



2 x 10- 4 Li 




) 



(i 



3 x 1Q- U cm, 



(20) 



where Lpi, R^ss are respectively the Planck length and the Schwarzschild radius of the earth, 
and are explicitly given by Lpi = 1.6 x 10 -33 cm, R&ss = 0.89cm. Then, it is of interest to 
notice that M* is approximately equivalent to the energy scale where new physics beyond 
Standard Model appears 



above which the coupling between the gauge field and the scalar field becomes strong. 

In this context, it is tempting to identify the it scalar field with a candidate of dark 
matters. Actually, this identification seems to be consistent with the results of OPERA and 
SN1987A at the same time by the following reasoning: In general, dark matters are expected 
to have a tendency of localizing near massive objects such as stars and the earth owing to 
gravitational interaction, compared to empty regions of outer space. The reason why neutrinos 
from SN1987A to the earth traveled at almost the same velocity as the true velocity of light 
in vacuum is that since there are not so much dark matters in outer space, the neutrinos 
from SN1987A propagated at the true velocity c without interacting with dark matters. On 
the other hand, since it is expected that there are sufficient dark matters on the earth, the 
interaction between the photons and dark matters reduces the speed of light on the earth to 
the smaller observed velocity c from the larger true velocity c. 

4 Discussion 

In this article, we have presented a resolution to one serious theoretical problem, that is, 
the bremsstrahlung effects of superluminal neutrions. Our resolution is on the kinematical 
grounds and thus strictly forbids the OPERA neutrinos to emit a pair of electron and positron 
via the bremsstrahlung effects. In this article, we have introduced the Galileon-type of scalar 
field to interact with the gauge field, but it would be possible to consider the other fields such 
as spin 2 new tensor field and spin 1/2 spinor field instead of the scalar field. 

It is then natural to ask ourselves if our resolution also provides a resolution to the other 
problems associated with the OPERA results. In particular, it is now known that the other 
challenging theoretical issue lies in pion decay process [3, 4, 5] where it was found that the 
decay of charged pion n + — > fi + + u^, which is nothing but the neutrino production process in 
the OPERA experiment, becomes kinematically forbidden for E v > 5GeV, which is obviously 
inconsistent with the OPERA results. It is worth stressing that our resolution also provides 
a solution to this problem since this decay process is automatically prohibitted whenever the 
OPERA neutrinos are subluminal. 



M* « ITeV, 



(21) 
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Finally, let us comment on the connection with related works. Moffat has presented a 
similar idea by taking account of bimetric relativity [17]. Li and Nanopoulos have recently 
advocated an idea that all the Standard Model particles might be subluminal due to back- 
ground effects on the basis of a string theory-inspired model [18]. Even if our present theory 
is very different from these works, it might be interesting to investigate the relation more 
closely in future. 

Acknowledgements 

This work (I.O.) is supported in part by the Grant-in- Aid for Scientific Research (C) No. 
22540287 from the Japan Ministry of Education, Culture, Sports, Science and Technology. 

References 

[1] OPERA Collaboration, T. Adams et al., arXiv: 1109.4897 [hep-ex]. 

[2] A. G. Cohen and S. L. Glashow, arXiv: 1109.6562 [hep-ph]. 

[3] L. Gonzalez- Mestres, arXiv: 1109.6630 [physics. gen-ph]. 

[4] X. -J. Bi, P. -Fei Yin, Z. -H Yu, Q. Yuan, arXiv: 1109.6667 [hep-ph]. 

[5] R. Cowsik, S. Nussinov, U. Sarkar, arXiv:1110.0241 [hep-ph]. 

[6] ICARUS Collaboration, M. Antonello et al, arXiv: 11 10.3763 [hep-ex]. 

[7] S. R. Coleman and S. L. Glashow, Phys. Rev. D59 (1999) 116008, arXiv:hep-th/9812418. 

[8] V. Gharibyan, Phys. Lett. B611 (2005) 231. 

[9] G. Dvali and A. Vikman, arXiv: 1109.5685 [hep-ph]. 

[10] L. Iorio, arXiv:1109.6249 [gr-qc]. 

[11] A. Kehagias, arXiv: 1109.6312 [hep-ph]. 

[12] P. Wang, H. Wu and H. Yang, arXiv: 11 10.0449 [hep-ph]. 

[13] E. N. Saridakis, arXiv:1110.0697 [gr-qc]. 

[14] I. Oda and H. Taira, arXiv:1110.0931 [hep-ph]. 

[15] A. Nicolis, R. Rattazzi and E. Trincherini, Phys. Rev. D79 (2009) 064036, 
arXiv:0811.2197 [hep-th]. 

[16] N. Nakanishi, Soryushiron Kenkyu 9 (2011) No. 3 (in Japanese). 



7 



[17] J. W. Moffat, arXiv:1110.1330 [hep-ph]. 

[18] T. Li and D. V. Nanopoulos, arXiv: 11 10.3451 [hep-ph]. 



8 



